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PERTURBATIONS OF 
NON SELF-ADJOINT STURM-LIOUVILLE PROBLEMS, 
■^ ■ WITH APPLICATIONS TO HARMONIC OSCILLATORS 

o 
o 

CN . LAURENCE NEDELEC ' 

Abstract. We study the behavior of the hmit of the spectrum of a non 
»vj . self-adjoint Stumi-Liouville operator with analytic potential as the semi- 

classical parameter h —> 0. We get a good description of the spectrum 
and limit spectrum near oo. We also study the action of one special 
perturbation of the operator (adding a Heaviside function), and prove 
that the limit spectrum is very unstable. As an illustration we describe 
t/J I the limit spectrum as ft ^ for P^ — —h^A + ix^ and the effect of this 

ri ■ perturbation. 



_. 1. Introduction 
CN ; 

^ , This paper is devoted to non self-adjoint Sturm-Liouville problems. We 

Qv I study the spectrum of the 1-dimensional, semiclassical Schrodinger operator 

■r^lj- ' on L^([— 1, 1]) with Dirichlet boundary condition, given by 

o 



X 



The potential y is a complex valued function on [—1,1], which extends 
holomorphically to some domain in C The boundary value at ±1 play no 
^ I special rule but are fixed to avoid more notation. 

S ■ The study of such operators is motivated by the Orr-Sommerfeld equa- 

^ I tion with linear profile 6, or by the non linear Zakharov-Shabat eigenvalue 

problem, cf. work of Miller jl.Sj . 

As an application we will focus on the case where 

(2) V{x) = ix^, 

and shall write P for the corresponding operator. Prom this one could also 
study the slightly more general case — ^^^^ + e'^x^, d G C using a change 
of variable. The spectrum of this operator on M was analyzed by Davies |2] , 
cf . also the recent work of Hitrik [9,^ . 

It is well known that the spectrum of a non self-adjoint operator is unsta- 
ble under perturbation of the operator. This motivates the introduction of 
the pseudo-spectrum, which has now been studied by many people, particu- 
larly Trefethen (who maintains the web archive http:web.comlab.ox.ac.uk/projects/pseudospectra) 



Supported in part by the FIM of ETHZ. 

1 



and Davies |S1)|1]; we note also the recent paper of Denker, Sjostrand and 
Zworski 0. 

We also consider the following perturbation of H: 

For P e (-1, 1) and d>0, let Hs,f3 

/Q^ Tjh 1,2 d ^ jr / \ f V{x) + i5, X > (3 

(3) Hs^P = -^d^ + ^^'/5(^)' ^^'/^Cx) = I ^y _ ^^^ ^ ^ ^ , 

with domain 

(4) Vs,^ = {uG L\[-l, 1]), u" e L\[-l, 1]), n(-l) = n(l) = 0, }. 

Notice that i/^^^ = H'^. 

Redparth ^^ has obtained similar results for arbitrary piecewise linear 
complex-valued potentials, a class which includes perturbations as above of 
the operator — /i^A + ix. 

We focus on the semi classical limit, and denote by lim/^ Sp{H^) the set of 
all values of E which arise as limits as /i — > of some sequence of eigenvalues 
of H^. 

We give a geometric criterion for E not to belong to lim/j Sp{H^): namely 
lim,, Sp{H^) C T/, T/ = C C T where Ts is the set of -E G C such that 
there exists a progressive path from —1 to 1 with respect to Vs^/s- ( We 
define this concept later) . This is done in section |3J see Theorems HJ This 
criterion is not sharp, and also not so easy to use. It can be proved using 
either ellipticity or exact WKB . 

As a second step we give an alternate criterion for E to belong to 
lim/i Sp(-ff^), when E is large enough, which is sharp and computable, see 
Theorem |21 . 

As an application we prove that the perturbation 6 changes the spectrum 
near oo quite drastically, see Theorem[7J The pseudo spectrum is defined to 
be the limit as 5 —> of the union of spectrum over all pertubations of H^ 
of size 6, and we show that this perturbation is sufficient to reach the entire 
pseudo spectrum of H^, see sectional We note also the work in progress of 
M. Hager that studies the effect of different pertubations on the spectrum. 

We use the geometric criterion for P to characterize the full set T'^, see 
Theorem llOL using special properties of the potential, and see that T^ forms 
a Y shape. After this article was completed, we learned that some of our 
results and methods are similar to those of Shkalikov [Hj , ^Sl ■ The algorithm 
proposed by Miller in 13^ to study the spectrum is the same as the one we 
use here, though he applies it to a different operator. 

We also obtain a description of the large spectrum of H which is near to 
lim/j Sp(P ) see Theorem El In the self-adjoint case; this can be found in 
the book of Marchenko JJ or the paper of Kappeler and Mohr JH]) and 
in some special non self-adjoint cases in work of Carlson, Threadgill and 
Shubin p. However, these results apply for potentials V with V' ^ Lj^c^ 
hence do not pertain to Vs^jj. 
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Even though it is not self-adjoint, the operator H^g has discrete spectrum. 
Indeed, Vs^/3 is compact relative to A, and from Weyl's theorem it follows 
that its essential spectrum is empty. In particular Hg^ has no residual 
spectrum. 



2. Generals tools 

In this section we are interested in the Dirichlet eigenvalues of the operator 
H^ on L2([-1, 1]) defined by 

f s j H^u = Eu, 

^^' \ n(-l) = n(l) = 0. 

Definition 1. A point xq in the complex plane is a turning point of order 
k for the operator H^ — E if V{x) — E vanishes to order k at x = xq. 

Definition 2. Let x be a point in the complex plane. The Stokes line 
starting from x is the set 

C{x) = {y e C, Re I {V{t) - E)^/"^ dt = 0}. 

J X 



This is independant of the branch of the square root. Denote by Sx y{E) = 

jy{v{t)-EYi^dt. 

Notice that the Stokes lines are integral curves of the Stokes vector field 



(6) s{x)=i{V{x)-EY/'^. 

The local structure of Stokes lines can be easily investigated (see e.g. J12j . 
0)- Away from turning points or singularities of V , the Stokes lines form 
a non singular analytic foliation. Indeed, if il C C does not contain any 
singularities of {V — E)^/'^^ then x h^ z{x,E) = J^ {Yif) — E)^/'^ dt is an 
analytic diffeomorphism. Notice also that from a simple turning point xq 
emanates three Stokes lines, each making an angle of 27r/3 with any other 
at Xq. 

Using ideas from quantum resonance theory, we shall consider as in [7] 
some distorted operators associated to H. Let 7 : [—1,1] 9 t — > C be a 
smooth simple path in the complex plane 7(— 1) = —1 and 7(1) = 1. Since 
V is analytic, we can define an operator H'^'^ on L^([— 1, 1]) as the restriction 
of H^ to 7([— 1, 1]). One computes that 

y\x) dx y[x) dx 

The correspondence between the Stokes geometry and spectral properties 
of H is given in the following result from 7~ . 

Proposition 1. The operator W^'^—E is elliptic if and only ifj is transverse 
to the Stokes lines. 
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Proof: The principal symbol of H'^''^ — E 



IS 



h'it^r) = j:^,{r' + W{t))\V{^{t)) - E)). 

We set {V{'~^{t)) — E)2 = re*^ and 7'(t) = e*^ (one can always suppose that 
|7'(t)| = 1). Then the path 7 is transversal to the Stokes lines if and only if 

det(7'(i),s(7(t)))/0, 

where s{x) is the Stokes field defined in Q. This condition is the same as 

Re7'(t)(F(7(i)-i?)'/V0, 

or finally 

cos(0 + e') ^ 0. 
Therefore, the transversality condition is equivalent to 9+9' ^ tt/2 mod [vr]. 
Since Im {-/' {tf h-y) = rHm2{9 + 9') and Re {-f'{tfh^) = t^ + r^ cos2{9 + 9'), 
the transversality condition is equivalent to the invertibility of K^, as stated 
in the proposition. ■ 

We deduce the following proposition, but this will be improved in The- 
orem n below, where we control the dependence on E and introduce the 
perturbation. 

Proposition 2. For h small enough, E ^ Sp{H^) if there exists 7 form —1 
to 1 which is transverse to the Stokes lines. 

Proof: General properties of distorted analytic operators imply that 
E G Sp{H^) is equivalent to £' G Sp(-ff'^''^). Let 7 be a path transverse to 
the Stokes lines. Then by Proposition ^ H"^' — E is elliptic and invertible 
for sufficiently small h (depending on E). ■ 

3. The semi-classical limit of the spectrum 

The section is devoted to the proof of the Theorem ^ and a direct appli- 
cation Theorem |2j The next sections H] and El are devoted to the proof of a 
converse of Theorem |2j 

We now introduce some notation. Denote by lim/j Sp{H ) the set of all Eq 
such that there exists a sequence hj — > with E^ . G Sp(i7''j) and E^ . ^ Eq. 
For a G C, we denote by zq{x) and zi rix) the action integrals 
(7) 

zo(x) = / {V{t) - Ef/^ dt, zi{x) = I (y{t) -i6- Ef/'^ dt, 

Zr{x)= f {V{t)+i5-Ef/^dt, 

where 7(— l,x) and Ji^ril^yx) are paths from —1 or /3 to a; on the Riemann 
surface associated to the determination of Sq of t ^^ (^(0 ~ EY'"^ and S;^,, 
associated to the determination of 1 1— > {V {t)+i5 — EY'"^ , S G C respectively. 



We denote by 7^ the set oi E £ C such that there exists 7i(— 1, /?) on which 
the function t i— > Re (2/(7/ (t))) is strictly monotone, and a path 7r(/?, 1) on 
which the function t >-^ Re (zr(7r(i))) is strictly monotone. We will call 
such a path 'progressive' . For 6 = we take a different definition: E £ Tq 
if and only if there exists a path 7 from —1 to 1 on which the function 
1 1— > Re {zoi'y{t)))) is strictly monotone. 

Finally, for E £ Ts, we denote by 

(8) d(^,r/,7) = inf inf |Re5^.z''''(x)| if (5 7^ 

l,r xe-y''^ 

(9) d(^,ro^7) = inf|Re9,z°(x)| 

xG'f 

where 7 is a progressive path associated to -E (7 = Ui_r7''^ if (5 / 0). 

Remark 1. Notice that 7^ is open. More precisely, if c is small enough 
there exists C > 0, such that if Eq belongs to Tg then, for E such that 
\E - Eo\ < cd(So,'7'/,7), we get E e Ts and d{E,Tf,-f) > Cd{EQ,T^',^). 
This follows from the relation 

VVs,p -E = y^Vs,i3 -Eo + 0(- ^"^" ^ 



which shows that a progressive path for Eq is also a progressive path for E. 

Theorem 1. Let 5 > be small enough and E eT^. Then E ^ Sp(i7^) as 
soon as h < d{E, T^, 7)^. (The reason for the exponent 7 will emerge at the 
end of the proof.) 

From Remark n and Theorem ^ we obtain the 

Theorem 2. If Eq e lim^ Sp(iJ^), then Eq E T/. 

Proof: Suppose that Eq ^ Tg and 7 is a progressive path associated to 
£^0- Let Eh be a sequence in the spectrum of H^ tending to Eq. By the 
remark above, E}^ belongs to Tg as soon as h is small enough, and we get 
d{Eh,T^,^) > Co- Theorem n shows that Eh is not in the spectrum of H^ 
as soon as h is small. This contradicts the hypothesis, so Eq G T^. ■ 

Proof of Theorem ^ We denote by (n'|_ ,uL) a basis of the space of 
solutions of the equation 

(10) (H^ - E)u = 0, 

in L^([— 1,/3]), and by (n^,n!l) a basis of the space of solutions of (jTU)) 
in L^([/3, 1]). A function / E L^([— 1,1]) is an eigenfunction of Hg with 
eigenvalue i? G C, if and only if /(I) = /(—I) = 0, and both / and dxf are 
continuous at x = (3. Writing 

(11) / = a^^u\ + aW_ on [-l,/3] / = a\u\ + aLuL on [/3, 1], 
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Then / is an eigenfunction of H^ with eigenvalue i? G C if and only if 



(12) 



/ uL{-i) 


nf,(-l) 








ui{(3) 


<{P) 


-uLiP) 


-^;(/3) 


d.u'.ifi) 


d,u^+iP) 


-d.u^-{(3) 


-5.n%(/3) 


\ 





U^_{1) 


ulil) 



f <\ 



J 



V«- / 



Therefore E belongs to Sp{H^), the spectrum of Hg, if and only if 
(13) detI{f3,E) = 0, 

where I{(3, E) is the matrix appearing here. 

Remark 2. For 5 = Q, the proof is slightly different, the corresponding 
matrix is 

n_(-l) n+(-l) 
u_(l) n+(l) 

We leave details to the reader and in the following, treat only the case 5^0. 



(14) 



I{f3,E) 



In the sequel, we shall compute this determinant for two particular bases 
(n'|_ , uL ) and [u^ ,u^-), namely for complex WKB solutions as in jH] , or j7j . 
Suppose E G T^. Then there is a path 7; = 7;(— 1,/3) transversal to the 
Stokes lines C{y), and we can suppose that 7i(— 1,/?) is of type +, that is 
t I— > z/(7/(t)) is strictly increasing. Indeed if it is not the case, the path on 
the other sheet of the Riemann surface S; associated to the definition of the 
square root with the same projection on C as 7;(— 1,/3) is then of type +. 
(We have zi[—l) < 0.) We define two independent complex WKB solutions 
w± of the equation {Hg — E)u = on the interval [—l,/3]. These have the 
form 

w^^:x^ {V{x) -id- E)-h^''^''^/''Wi{zi{x)) 
(15) 

= (a,.z,(x))^2(x)e±^'(")/'^H^i(z,(x)). 

Similarly, there exists a path 7^ = 7r(/?, 1) transverse to the Stokes lines 
of type + , (We have Zr{l) > and we have again two independent complex 
WKB solutions of the equation H^w± = on the interval [/3, 1] : 



(16) 



W4 



{V{x) + i6- E)-h^'-^^^/''W^{zr{x)) 

= {d,Zr{x)rh^'^'-''y^WUZrix)). 



rLr 



The functions Wl are convergent series of the form 



(17) 



W'i' 



n=0 



±! 



rUr 



where the W^j^ are determined by the recurrence relations: 



(18) 



{dz ± l)W, 



l,r 
2n+l,± 



l,r-iitrl,r 



-W'^'W. 



2n,± 



d^W^ 



l,r 
2n,± 



hrur^^-r 



-H'^^Wr 



2n-l,± 



with initial data 



( <± = 1, 



(19) 



wl_,{zi{-i)) = 0, vf;_(z.(i)) = 0, n > 1, 

Wlizim = 0, WlAz^m = 0, n > 1, 



We have set here 



(20) 



H^'^'iziA^)) 



1 d^V{x) 



4 (9,^i,r(x))3 
These equations can be written in integral form as 



(21) 



where 



(22) 






l'£{v){z) 



ll,r(z) 



o±'ii'---)/hH^^-^\u)v{u)du 



.r^^{v){z) 



ll,r{z) 



H''''{u)v{u)du 



and 7«,r(-2) is the image by x h^ zi^^i^) of 7; ,, from the initial point — 1 to z 
for I and from P to z for r. 

Now we have the following estimates 

(23) 

sup \I^/{v){z)\ < sup \v{z)\ sup {l^''''(^)l | '!f':^^,l!,| }x 



^^7; . 



2(n — z) -, 



|Re3a;Zi,r(a:^) 



J-y f ZiU — Z) \ 

5R^(^ir^)Re(dz)| 



7l,. 



< /i sup |t;(z)| sup |ff''''(2)||a^z;^^(x)| — - 

< C{V)h sup |z;(z)| ^ 



d(^,'r/,7)3 



(24) 



sup I J '''(f)(2;)| < sup \v{z)\ sup |i:/' '^(z)|(92:Zi^r(a^)| 

< sup \v{z)\ sup |F''''(z)||a^z/^r(a;)| 

ze7i,r ze7,_r,2=z'.''(a:) 



< C{V) sup |f(z)| 



1 



d(.E,r,-,^y 



Similarly we obtain 

(25) 



rl.r. 



1 



sup \lf{v){z)\<h sup \v{z)\ sup \H '''iz)\\d^zi^rix)\-, , 

We denote by [xj the integer part of x. Hence on a progressive path 7 we 
have 



l<',±|oo 



< 



ci(i^,r/,7)3-' 



which means that 



c ^,\H )j 



>7) 



(26) 






d(B,r/,7)'s)- 



The function w^ form a basis of solutions of the equation Hsw = 0, and 
similarly w!^ is a basis of solutions of the equation Hgw = 0. A computation 
gives 



Lr 



"'"'^'-(x) = ±a,z,,,(x)5e±^i j;(-l)"<'^(z;,,(x)), 



dx 



n=0 



and so 



Lr 



dwX. 1 I ^i,r(^) 1 /l 



))■ 



Now compute the determinant; we shall write O = 0( ,.^^e — vf) for sym- 
plicity. 

(27) 

/ wl{-l) u;V(-l) \ 

u;L(/3) <(/3) -u;:(/3) -<(/?)) 

' ti;''_(/3) u;''+(/3) -«;'"_ (/?) -w^'^P) 

V u;1(l) ti;';(l) / 

1 -i -i -i -i ZjC-l) ^r(l) 

a^j^r '(l)5x-2r '(/3)9^Z; 2(/3)a^2;^ 2(-l)e h e h X 



/l 














1 + 

1 1 + -1 + -1 

-a,zK/3)(l + 0) 5,^K/3)(l + 0) 5,^r(/3)(l + 0) -a,Zr(/?)(l + 0) 











1 + 



So the sign of the determinant is given by the sign of 



f I o\ 

(det 



1 -1 

d.ziip) d^zriP) 

\ 1 / 



^^(i(E,r/,7)6^^ 



(28) 



= d^Zi H-l)d^Zr Hf3)d^Zi ~'{(3)d^Zr ^(l)x 

az,(/3) + d^zm + 0(-— A^— )) 

= {V{-1) -i5- E)-^V{f3) -i6- E)'^ 
X (V{p) +i6- Ey^V{l) +i6- E)'^ 



h 



X {^V{f3)+i6-E+ ^Vil3)-i6-E + 0{ ^^^^^^^^^J ) 



( we have Re {^/V(J3)+i6^^ + y^V(J3y^i6^^) > 2d(^,T/,7).) 
The determinant is non zero as soon as ,,i-,^c — \Tr is smah and als 
compared to d{E,T^,^) (and this explain the exponent 7). 



4. SPECTRUM FOR LARGE E AND 6 = 

In this section we prove three results: the first, Theorem |21 gives the 
condition for £"0 to belong to lini/i Sp{H ). The second, Theorem|3J describes 
lim/i Sp{H ) as a curve and gives its asymptotics near 00. The last. Theorem 
El describes the eigenvalues of H which are near to lim/iSp(iJ ). In the 
following we let denote by y a primitive of the potential V. In this section, 
we assume that V verifies the following hypothesis: (HI) For any fi C C, 
if V~^{Q) bounded then Q is relatively compact. 

(H2) If E is large enough, but with small imaginary part, then 1 and 
— 1 belong to the same Stokes region. We recall that this means that one 
can find a path going from —1 to 1 that does not intersect the Stokes lines 
issuing from the turning points. Notice that, under the assumption (H2), 

E £ Tq is equivalent to Re (/_-^ \/V — E dx) = 0. 

Theorem 3. For Eq large with small enough imaginary part, then 
Re ij^^ VV - Eq dx) = if and only if Eq £ lim,, Sp(i/^). 

Moreover there exists C small enough and Eh in the spectrum of Hq which 
satisfies |\/S^ — \/^| < Ch. 

Proof of Theorem \^ The reverse implication follows directly from Theo- 
rem |2l so we prove the direct implication. 
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There exists two WKB solutions w^ of Hq — E = of the form 

(29) w^{x) =d^zo{xyh^"-'^''^/''W±{zoix)), 

W±=Y,Wn,± H^o,±(-l) = l, VFn,±(-l) = Oforn>0. 
i=o 

The complex number E belongs to the spectrum of Hq if and only if 



^0(1) 



W+{1) 



which is equivalent to the existence of /c G Z such that 

/I J-i 1^ 

^V{x)-Edx - -\n{W+{l)) + -ln(VF_(l)) - ikkir = 0. 
-1 2 2 

ate 

f{E)= I yJV{x)-Edx-lm I ^/V{x)-Eodx, 

Im f\y/V{x)-Eodx 



Let us denote 



I/O ux, 
and 

ko=[ 



and also 



/I i-, i-i 

y/V{x)-Edx - -ln(l^+(l)) + -ln(H^_(l)) - ihkoTT. 

We have f{Eo) = 0. We want to apply Rouche's theorem, but to do so we 
must give upper and lower bounds for Vl/'-[_(l) and an estimate on /'{Eq) = 

By the hypothesis on the geometry of the Stokes lines near —1 and 1 there 
exists a path 7 which links —1 to 1 on which Re J_^ y^V{x) — Eq dx = for 
all t G 7. Moreover ^^(7) is bounded uniformly with respect to Eq, and by 
assumption (HI), 7 is also bounded uniformly with respect to Eq. 

We have 

' ^(1 + 0(2^)), 



therefore 



^Vix) - Eq i^o' Eq 

\f{E)-f{EQ)\>-^\{E-EQ)\, 
2vl^o| 

sup |V^'(a;)| 

sup\H{z{x))\ <Ci 3 — , 



and 



7 |£^o|2 



fb' fp — fp 

sup |Re / yJV{t) - E dt\ < C2 -° , 

i,fe'67 Jb V 1-^0 1 



b,fe'67 
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for any E and ag such that \E — Eq\ < o(Eo) and \Eq\ > Cq . Using the 
expression of Wn,+ {z) as a Volterra integral, we obtain the estimates 

2 r'^' 1 

\Wn±(z)\ <exp|- sup iRe / ^/V(t) - E dt\(n + 2)| sup \Hr — . 

hb,b'£-y Jb ' -y n\ 

This gives 

C2 \E — Eq\ -, C" 



\Wn,±iz)\ <exp{-^^=fn + 2} 



for \E — Eq\ < ^o(-E'o) and |£^o| > Cq • So we obtain 

W±{l) = l + -—Te^v{-^ ^^}exp{— ^e >^ ]. 

This gives the estimate 

\W^(l) _ i| < _^eC2C3gCi|i=;orie^2C3 
|^o|5 



for Eq > Cq and \E — Eq\ < Ca/iy^l-Eol with h small enough. We also have 
(30) |lnVF+(l)| + |lnVF_(l)| < V 

\En\^ 



for £"0 > C{Cs), and |-E — £'o| < Cs/iy^jE'ol with h enough small, and for 
the functions / and g we get 

\f{E) - g{E)\ <h + ^(|lnl^+(l)| + |lnH^-(l)|). 



Thus if \E — Eq\ = C^hy^lEol, for some large constant C3 > 0, we have 

(31) \f{E) - g{E)\ < (1 + C,\Eo\-l)h, \fiE)\ > ^h. 

Finally, if Eq is large enough, all the assumptions of Rouche's theorem are 
fulfilled, and we get the result. ■ 

Theorem 4. For any fixed a large enough, the equation 

Re5_i,i(£)|^_^^^ = Re / ^/V{x) - (a + ib) dx = 

has a unique solution 6(a). Moreover 6(a) = ^^ra (^(1) — Y{—1)) + 0{-). 

Proof: First we prove that if 6 is such that Re J_^ ^yV{x) — a — ibdx = 0, 
and 6 = o(a) then 

(32) b = i-lm{Y{l)-Y{-l)) + 0{a-^). 

Indeed if we denote by -E = a + i6 

ry 

ip{E,a,y) = / ^/V{x) - Edx, 

J a 
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then, uniformly for y in a compact set, 

^{E, a, y) = iVEiy - a) - i-^{W{y) - W{a)) + 0{E-i). 

ly E 

Since Re5_i_i(-E') = 0, we have 

Re {i2^ - ^^(^(1) - ^(-1))) = 0{E~'i). 
Writing yE = c + ib, then b = o(c) and 



1 „, 6 

-26 + 

or equivalently 



26 + j-lm {¥{!) - y(-l)) + O(^) = 0(c^3) 



6 = llm(y(l)-y(-l))+o(i). 
4c c 

This gives 6 = ilm (y(l)-y(-l))+o(l), and 6 = 0(1). Therefore 6 = 0(i) 
and recychng through this argument, 

6 = lim(y(i)-y(-i)) + o(i), 

which gives (|32|). 

The existence part of the theorem is straightforward: set 

ip{c, 6) := Re / \/cV{x) -l-ibdx = cRe / Vy(x)-S(ix 

where E = i(l + i6). Since db^p{0,0) = -2 and ((3(0,0) = 0, the imphcit 
function theorem apphes. 

Finahy, the uniqueness fohows from the fact that the map (p : b ^^ 
Re J_^ Y^y(x) — a — ibdx is injective for a large enough. Indeed 



\^/V{x)-a-ib- y^V{x)-a-ib'\ < C^y\b-b'\, 
for a suitable branch of the square root. ■. 

Remark 3. We can compute more terms of the asymptotic expansion than 
in If32\) . Indeed, from 



(33) 



1 r^ 

ip{E, -1, 1) = i2y/E - i^^{Y{l) - Y{-1)) + iE-2 / V'^{x) d. 

-iE-^ / V^{x)dx + 0{E-2] 



we get 

b{a) = iim (y (1) - y(-i)) + Aim (1'^ V') + o(l). 

Looking more carefully at the quantization rules, we obtain the asymp- 
totics of the eigenvalues. 
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Theorem 5. Let Eq be a solution of Re ( J_^ \/V — Eq dx) = with \Eq\ 
large and Im {Eq) small. If E^ is in the spectrum of Hq and satisfies \^/Eh — 
^/EqI < Ch, provided C is small enough then E^ satisfies: 

^TThk^, , (y(l)-y(-l)) (y(i)-y(-i))2 i 



2 ' 2 {2hkTTY \hkY'- 

for some /c G N. 

Remark 4. For a real potential V hounded on [—1, 1] and for any Eq ^"K 
large we get Re (/_-|^ yJV — Eq dx) = 0. So Theorem\^ can be read as: If Eh 
is large and satisfies Eh G Sp{Hq) then Eh satisfies: 

^_^^hk^, ^ (y(l)-y(-l)) ^ (y(i)-y(-i))2 ^ ^^ i ^ 



2 ' 2 {Ikk-n)^ \hk)^'' 

for some /c G N. 

Proof : 

/•I 1 3 

/ VV-Edx = i2VE - i^^{Y{l) - Y{-1)) + 0{E-^). 

Then E G Sp{Hg) if and only if there exists fc G Z such that 

/I P, Pi 

^V{x)-Edx - -ln(W+(l)) + -ln(Vr_(l)) - ihkir = 
I 2 2 

For Eq big enough, \E — Eq\ < C^h^/EQ, C3 is as in (pTjl and h smah, we 
get 

|lnVr+(l)| + |lnVF_(l)| < Cl^ol"^ 

So -E G Sp(/7^) if and only if there exist fc G Z such that 

z2/B - i-^{Y{l) - Y{-1)) - ihkiT = 0{E-^) 
1y E 



In particular taking \/ E = c + ib the real and imaginary parts give 
2c-hk7r = 0{c'^) c=-hkTr + 0{{kh)~^) 

2b = ^Im (y(l) - y(-l)) + 0(c-3) 
2c 
Using this again gives 

c = \hkiT + ^Re (y(i) - y(-i)) + 0{{kh)-^) 



Remark 5. One can also treat the case P on all o/M. Then a condition 
for E to belong to lim/j Sp(ff ) is that there exists a progressive path joining 
two turning points. 
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5. Large spectrum of a perturbation of H 

In this section, we prove two theorems: Theorem |H1 gives the condi- 
tion for Eq to belong to hm/j Sp{H^) and Theorem [7| gives a description 
of hm/i Sp{Hg) as the union of two curves, the asymptotics near cxd of which 
are made exphcit. 

We assume that V verifies the hypothesis (HI) and in addition : 
(H3) If E is large enough, with small imaginary part then 1,-1 and (3 
belong to the same Stokes region. 

We assume /? 7^ ±1. Recall that (H3) means that one can find paths 
going from — 1 to /5 and from 1 to /? which do not intersect the Stokes lines 
issuing from the turning points. With the hypothesis (H3), the condition 

E £ T^ is then equivalent to 

either Re (/f^ VV - i5 - E dx) = or Re (/j yjV + i5-Edx) = 0. 

Theorem 6. For Eq large enough, with small enough imaginary part, then 
Eq satisfies Re (J'_-^ y/V — E^dx) = or Re ( L ^/V — Eq dx) = if and 
only if Eq £\imhSp{H^) 

Theorem 7. For any fixed a large enough, the equation 

Re (5ii JE))\E=a+^b ■■= Re / yJV{x)±id-a-ihdx = 
has a unique solution b[a), and this solution satisfies 

h{a) = ij^^^ (^(/3) - ^(±1)) ^i^ + 0{\). 

Proof : The proof is the same as in Theorem ^ where the potential is 
V ^i6 in each side and we just need to change the interval for the integral. 
■ 

Proof of Theorem^ The reverse implication is already proved in Theorem 
1^ Let us prove the direct implication as in the proof of Theorem |3J There 
exist four WKB solutions w^ of H^ — E = Q on the interval [— 1, /?] or [/3, 1] 



(34) 

with initial data 

(35) 



u-^ = (a^Zi(x))-^e±^'(^)/'^iy4(zKx)), 
w;!t = (5^z^(x))-ie±^'-(^)/'*WX(2r(x)), 



<± = 1, 



Wl^+{zim = 0, VF;„(zr(/3)) =0, n > 1, 
Wl_{zim = 0, Wl^+izrm =0, n > 1. 
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Recall that E G Sp{H^^') if and only if 



(36) 



= det 



/ wLi-l) w^_^i-l) 

w^'_{l3) w%{(3) -w'''_{l5) 
\ wLil) 



\ 

<(1) J 



-d^Zr ^{l)d^Zi ^(-l)x 



I e h Wi(z;(-1)) e h wi(z;(-l)) 



det 



Vo 



1 




1 



d:czf{l3) 

Ml) 



d^z?(f3) 
d^z?{P) 
d.zf(p) 

r-(l) 



"W!:(zr(l)) e~?^W:^(zr(l)) / 



To make the computation, simplify the notation by setting 



d:,zf{(5) 



zr(l) 



Hi-^) 



^ y = e h X = e h 



d^Zrm3) 

then (|5S)) becomes 

(37) 

/ x~^WL{zi{-l)) xW[{zi{-l)) 

det 



\ 



which gives 
We remark that 



1 
1 




0, 



yWLizril)) y-^Wlizr{l)) ) 



^ WL{Zr{l))l-t' .,l-f 2 WL{Zr{l)) 



-1 



So the condition for E to be in the spectrum of H^ is 



(38) 



^H^ w[{zi{-i)) _ {d^zrm-d^zim ^ 

Wl{zi{-1)) {d^Zr{[3) + d,Zim' 
_^._^ WL{Zr{l)) {d,Zr{(3) - d^ZliP)) 

Wl{Zr{l)) {d,Zr{l3) + d^Zl{f3)) 

. {d^ZriP) - d,zi{f5)f . 

"^ {d^Zr{(3) + d,Zi{f3)Y'- 
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) 



We write this condition by taking logs as 

2Mzi)T^i(2K-l)) {d.Zr{P)-d^Zi{p)). 



ln( 



e 



WL{zi{-l)) {d,Zril3)+d,Zi{P))' 



f . . .2^rm WL{Zr{l)) id,Zr{(3)-d,Zii(3)) . 

-in(i- ;f-;^|-f-;fii: )=2..., kez. 

^ {d^zr{p)+d^zi{p)f' 
Stepl First assume that 

Re ( / y^V -i6- Eq dx) = and Re ( / ^JV + i5- Eq dx) / 0. 

(The case Re (/^^ y/V - i5 - Eq dx) / and Re {j^ y/V + iS - Eq dx) = 
is treated the same way. ) 
As (jSni), we get 

(40) \lnWi{zi{-l))\ + \lnWL{zi{-l))\ < Cil^ol"^ 



for \E - Eo\ < C^hy^lEol and Eq > C big, h small. As in ^^, we get 

(41) |lnP^;(z^(l))| + \lnWL{zr{l))\ < Ch 
with C uniform in Eq, for \E — Eq\ < Cs/i y|^Eo| and 

/ . N , n ^r(l) , C 

(42) |e~ ~^| < e""h. 
Rewrite the previous condition as 

■ Wiizii-l))^ 



gi^k{E) = 2zii-l) + h\n{ 
+ Mn(l - e 



VFi.(zK-l))' 



(43) +Mn(§f4§^fMM) 

(9^Zr(/3) +5a;zK/3)) 



+ Mn(l + e~ 



Ty|(zi(-l))(9.z.(/3) + 5,zK/3)) 

.Jim. 
xzmy 

2M1) vFi(z..(i)) (a.z.(/5) + a.zK/?))^ 



vF;(z,(i))(9,.z,(/5)-a,zK/?))' 



Setting 



and 



,, , id^ZrW) - d^ziip))^ , 



fi{E) = 2zi{-l) + 2lm{ ^/V{x)-i6-Eodx) 



P 



[lm{l ^/V{x)-i6-Eodx)h~^\, 
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we have fi{Eo) = 0. Using ^, dHJ) and ^ gives 

\9i,koiE) - h{E)\ < h{l + C^5 + ChV5'^ + Cil^ol"^)- 
Now recall that 

\f,{E)-h{Eo)\>^=\{E-Eo)\. 
y\Eo\ 

We can apply the Rouche's Theorem to prove the existence of an eigenvalue 
of H'^ for each /i at a distance C3 {5) \/\EQ\h of £"0 for Eq large. Step2 Assume 

now that 

Re(/ yJV -i5- Eq dx) = Re ( / \/F + z5 - £"0 dx) = 0. 

As in (jSOl), 

|ln<(zK-l))| + |lnH^i(zK-l))| + |lnH^;(z.(l))| + 

\\nWL{zr{l))\<C^\Eo\-'^. 



For Eq > C, C sufficiently large, and l^ — -EqI ^ C's^vJ^ol with h small. 
We write dSH) as 

g,,,{E) = 2z,(-l) - 2zr{l) + Mn(^±M^) 

W^i(2;i(-1)) 

. .2^ wL{zi{-i)) {d^zm - d^zim 
(45) ''^ ' wi{zi{-i)) {d.zm + d.zmy 

and write 



^1^(1 - t^!ta^ , ^■,^^^2 - 2^^^^ = 0> 



/2(£) = 2zK-l) - 2zr{l) - 2Im( / ^V{x) - iS - Eq dx) 
(46) ^;, 

Let 

ko = Llm( /" V^(^) -iS-Eo dx - [ ^/Vix) + i6 - Eodx)h~^\. 

Jl3 Jp 

Then 

\f2{E) - f2{Eo)\ > ^Ae^\~^\{E - Eo)\, 



and for \E - Eq] < C^hy^l 



\h{E) - 52,fco(^)l < KC9\Eo\—^ + CsVd + 1). 
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On the set where \E — Eo\ = C^hy^jE^, we have for C3 = 03(6) and £"0 
large 

\{f2-92,ko){E)\<\f2{E)\. 

All the hypotheses of Rouche's theorem are now satisfied and we conclude 
as before. ■ 

6. Application to the harmonic oscillator 

In this section we are interested in the Dirichlet eigenvalues of the operator 
P'^onL2([-l,l]), 

(47) P'^ = -h'— + ix^ I ^''^ = ^^' 

^ '' ^ dx^^ ' \ u{-l)=u{l)={). 

We obtain here only two results. Theorem llfll describes the shape of the 
set lim/i Sp(P ), cf. figure ITU Theorem ll2l shows how the spectrum changes 
near 00 when P is change by a specific perturbation of size 5, as illustrated 
in figure El • 

For any E ^ C* there are two turning points, a±{E) = ^{—iEY'"^, which 
are both simple i.e. of order 1. 

Denote by Sx,y{E) = J^iit'^ — E)^/'^ dt. The function Sq^x is even, so to 
simplify the computation one can use that S^x,x = 25o,a;- Writing 



{x,y,E)= r {it" - Efl^ dt, 

J X 

'^x,y,E)= r {u^ - E)'/^ du, 

J X 



and changing coordinate t = e s u, we get 

z{x,y,E) = ±z"-"-^[xe>i ,ye>i ,Ee 4). 

So the Stokes lines of —h?IS. + ix^ can be deduce from those of the harmonic 
oscillator —h?/\ + x^ by a rotation by — | about the origin, jj^l- Using the 
geometry of the Stokes lines we find for some E a progressive path from — 1 
to 1. Combine with previous theorems, we obtain a partial description of 

Theorem 8. We have 

• Suppose E £ C* is such that ReSa^,a+{E) 7^ 0, ReS'_i^i(-E) 7^ 0, 
ReSa+,i{E) / 0, ReSa_,i{E) ^ 0, then E eT 

• If E £ T, then for any h small enough, E is not an eigenvalue for 

ph 

• IfE£ hmh Sp{P'') then E £ T'' 

Proof: Using Theorem ^ and Theorem |2l it is enough to find a path 7 
tranversal to the Stokes lines from —1 to 1. We have ReSa_,a+{E) / 0, 
so the Stokes line issuing from a+ do not intersect the Stokes lines issuing 
form Q_. Therefore the complex plane is divided into exactly five region 
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delimited by the Stokes lines issuing from the turning points. In Figure 
we have drawn the configuration of the Stokes lines up to an analytic 
diff eomor phism . 





Figure 1. Non degenerate Stokes lines 

We have Re Si^a±iE) ^ 0, so that in particular 1, (and by symmetry 
are not on any boundary of the Stokes regions. 




Figure 2. Progressive path issuing from 1 

If 1 and —1 belong to different regions, then there always exists a path 7 
tranversal to the Stokes lines (see the figure |2) ; if on the other hand 1 and 
— 1 belong to the same region, but we also have ReS-i^i{E) ^ 0, then they 
do not belong to the same Stokes lines and there still exists such a path 7. 
■ 

If ReSa_,a+{E) = 0, the Stokes configuration is as in Figure El 
We now picture the evolution of the Stokes lines as E moves in the complex 
plane. 
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Figure 3. Degenerate Stokes lines 




Figure 4. Stokes line for different E 

In fact Theorem |H1 can be improved. Choose the segment [a_ , a+] as a 
cut for X 1-^ {y{x) — E)^'"^. The pictures [3 and IHl show the branch of the 
square root for different values of E. 

Denote by 

r_i,i = {Ee C;ReS-i,i{E) = 0, Re5,_,«^(S) / 0} 

The next theorem is proved by listing all the situations where there do not 
exist progressive paths from —1 to 1 using the geometry of the Stokes lines. 
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-Re 



ill/8 




i Jt/8~ -- - , 



Figure 5. branch for E = e*4, the cut is the dashed hne 
and also a Stokes hne 




it/s i p;2 



Figure 6. branch for E = e*4e*^,/3 > smaU. The cut is 
the dashed hne 



Theorem 9. We have 

• Suppose E £ C* belongs to 

ra_,a+ ={Eg C;ReSc,_,^4E) = 0,ReSa+,iiE) < 0}, 

then E gT. 

• Suppose E £ C* belong to 

ra+,1 = {EG C;Re5„+,i(^) = 0;Im5„+,i(E) < 0}, 

then E gT. 

• Suppose E £ C* belong to 

ra_,i = {E£ C;ReSa_,iiE) = 0;lmSa_,iiE) > 0}, 

then E gT. 

• E ^ Ta^^i U Ta_^i U Ta_^a+ U F-i^i then there is no progressive path 
form —1 to 1. 

Remark 6. If a progressive path enters a Stokes region crossing a Stokes line 
issuing from a turning point a, then it cannot leave this region by crossing 
any others Stokes lines issuing from the same turning point a. 
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The proof consists in listing the cases where progressive paths do not exist. 
Recall first that if 1 and —1 belong to the same region then there exists a 
progressive path if and only if ReS-i^i{E) / 0. We get r_i i = {E £ 
C; 1 and — 1 belong to the same region Re5_i^i(£') = 0}. Now assuming 1 
and —1 belong to difi^erent regions, we get two different figures for the Stokes 
lines either ReSa_,a+{E) = (FigurfE} or ReSa_,a+iE) / (Figure HI). In 
the first case, using Remark El we see that we cannot find a path if and only 
if 1 belongs to the hatched region ( and —1 by symmetry to the opposite 
one) (Figure Ej) . 





Figure 7. The position of 1 is in one of hatched regions 



The condition ReSa^,a+iE) = 0, KeSa+,i{E) > corresponds to this 
situation, i.e. Figure El shows the sign of the quantity ReSa^,x{e^^^) de- 
pends on the position of x for A G M. 




Figure 8. Sign of the quantity j; ^ Re Sa+,x{e^''* , 



Now assuming 1 and —1 belong to different regions and ReSa_,a+iE) ^ 
(Figure El). Using Remark El we see that we cannot find a path if and only if 
1 belong the dotted curves ( and —1 by symetry to the opposite one ) (Figure 
ini). The condition ReSa+,i{E) = 0, IicaSa+,iiE) > correspond to Figure 

El 

The condition KeSa+,i{E) = 0, lm.Sa^,i{E) > corresponds to this 

case. i.e. Figure ITOl shows the sign of the quantity Im Sa±,x{s^^^ g^^) with 
/3 > and small depending on x near to a± . 
Now we want to describe the set T^. 
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Figure 9. The position of 1 is in one of the dotted Une, that 
are Stokes hnes 




Figure 10. Sign of the quantity ImS'Q,j^,j:(e* 4 Ae*^) 



Theorem 10. We have 

(1) The set T^ is the union of three curves Fc^^i, Ta_^a+, r_i,i with 
only one of inhnite length r_i^i. 

(2) The three curves meet at a common point Aqc^ with Aq G K^. 

(3) The curve r_i^i goes from Xoe~ to M + |. 

(4) Ta_,a+ = {AeX; A G M, < A < Ao}. 

i7V_ 

(5) The curve Ta^^i goes from i to Xqc * . 
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Proof: (1) : 

Let E = X^e^l^ with /? G M and A G R. We fix /3 and increase A. We get 
z{x, y, A^e*^) = z{j, j, e'^). The shape of the Stokes lines remains invariant 
up to dilation by A. Let (3 be such that ReS'a_,Q,^(e*^) 7^ then under 
symmetry we get Re5o,Q+(e''') 7^ , and ReS'a_.o(e*^) / 0. So if /? is 
such that there is no Stokes line issuing from the turning points and going 
through 0, then there exists a neighborhood of with no point of any Stokes 
lines issuing from the turning points. For A big enough, we get that —1 and 
1 belong to this neighborhood of 0. Then E will belong to T'^ if and only if 
Re5_i,i(S) /O. 

If (3 is such that Re5a_,„^(e*^) = then /3 = ±f and Figure [HI and the 
fact that the spectrum is included in the set {E G C; £' = M^ + ix'^,x G 
[—1,1]} (i.e. the pseudo spectrum or the values of the principal symbol) 
shows that E G T'^, E = Ae*4^ A G M, implies E bounded. So we have that 
if E is big enough and belongs to T"^, then E G r_i^i. 

We choose the determination of a± so that the set 

ra_,i = {^ G C; ReSa^,i{E) = 0; lmSa_,i{E) > 0} 

is empty see Figure^ 

The second statement is proved by the relation 

So,a4E) + Sa+,l{E) + Si,o{E) = 0, 

and illustrated in Figure ITTl 





Figure 11. Stokes line for E near Aqc 4 , the arrows locate 
the positions of —1 and 1 

Point 3 is proved in Theorem 1111 below. 
Point 4 is deduced from 

ITT 

zix,y,E)=± / ^^ {u^ - Ee-'^)^/Uu 

Jxe"S" 

So for E = {\e~ , A G M} we obtain 

Re 5o,a+ {E) = Re ± / {u^ - A)^/^ du = 0. 
Jo 



'0 
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The existence of Aq is obvious in Figure El 





Figure 12. Stokes line for E = \e ^ 
Point 5 is obvious from Figure [T3l 





Figure 13. Evolution of the Stokes line for E G Fa+.i, the 
arrows locate the position of —1 and 1 



The hypothesis (HI) is obviously satisfied by the potential V{x) = ix^ . 
The hypothesis (H2) (H3) are satisfied by the potential ix^ ( and by any 

potential e'^x'^ if ea ^ M ) but not by the potential x^. 

For large energies, we summarize the result of Proposition |2, Theorem 13] 
and Theorem m in 



Theorem 11. We have 

(1) Large values of lim/j SpP^ are close to the curve r_i^i. 
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(2) The curve r_i^i = {E £ C,ReS'_i_i(£') = 0} tends to infimty, and 
is asymptotic to the line M + i/3. 

One could compute the set T^ with Matlab or Mathematica. In Figure 
1141 we have drawn the set T^ using Theorems IIUI and 1111 



1 










/ \ ■'« 1 1/3 






/ ♦ 














Figure 14. Limit spectrum of P'^ 

For large energies and perturbation , we summarize the results of Theorem 
\^and Theorem\^in 

Theorem 12. We have 

(1) Large E G hm^ SpPj^ are near to the union of the two curves r_i ^, 

(2) The curve T_i^p = {E £ C,ReS'_i^^(-E) = 0} tends to infinity, and 
is asymptotic to M + i|(/?^ - /3 + 1) - i6. 

(3) The curve r/3^1 = {E G C,Re5'/3^i(-E) = 0} goes to infinity, and is 
asymptotic to M + ii(/3^ + (3 + l) + i5. 

We remark that the two curves r_i^^ F^^i are distinct for /? 7^ 0. 
Fisure [731 represents the two Theorems. \11\ and \12l 

7. PSEUDO-SPECTRUM 

We have defined here the pseudo-spectrum of Hq as the set 



{e + Vix),x£ (-l,l),eGM,Im(eF'(x))/0}. 
We remark that for any z is this set, we could contruct a function u in the 



domain of Hq such that 

\\{H^-z)u^\\ = 0{h^)\\u% 
Let suppose now that we put two jumps in the potential, one at j3 and 
one at /?'. i.e Vppi = V + i5\Hp + i62Hpi. 
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2 
P+P 



2 
P -P 



Figure 15. Spectrum of P in plain lines and spectrum of 
Pj for 6 small near to oo as dotted lines. 



Assume that V verifies the hypothesis (HI) and : 

(H4) If E is large enough, with small enough imaginary part then 1,-1, 
P and (3' belong to the same Stokes region. 

Then we obtain just as in Theorems\^and\^the next two theorems 

Theorem 13. For Eq large, the following two conditions are equivalent 
(1) Re il"^ ^Vp^>-Eodx = 0) or Re {f^, ^V^p.-E^dx = 0) or 

0). 



Re (// v^ 



Pfji — Eq dx 
h 



(2) EoelimhSpH^. 
Theorem 14. For a large enough, the equation 



^eSj3^pi{E)\E=a+ib 



P' 



Re / yJV{x) + i5i 



has a unique solution b{a), the solution satisfies 

1 



b{a)-- 

We remark that 
lim 



^^_^lm{Y{P')-Ym-i6, + 0{-). 



ibdx = 



.1. 



1 



/3'-/3 P'-(3 



lm{YiP')-Y{(3))=lm{Vm 



This means, letting f3' -^ [3, and 5 — > any values of the form {^^ + 
ilm.{V{x));^ € M, |^| >> l,x £ [—1, 1]} is in the spectrum of this kind of 
perturbation of H^ . Thus we can reach all the large values of the pseudo- 
spectrum with this special type of perturbation. 
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